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Abstract 

The magnetohydrodynamic (MHD) equations of incompressible viscous fluids with finite 
electrical conductivity describe the motion of viscous electrically conducting fluids in a mag- 
netic field. In this paper, we find twelve families of solutions of these equations by Xu's 
asymmetric and moving frame methods. A family of singular solutions may reflect basic 
characteristics of vortices. The other solutions are globally analytic with respect to the spa- 
cial variables. In particular, Bernoulli equation and Wronskian determinants play important 
roles in our approaches. Our solutions may also help engineers to develop more effective 
algorithms to find physical numeric solutions to practical models. 

Keywords: MHD equations; asymmetric condition; moving frame; Bernoulli equation. 
AMS Subject Classification (2000): 35Q51, 35C10, 35C15. 

1 Introduction 

In physics, magnetohydrodynamic (MHD) studies the motion of electrically conducting fluids 
in a magnetic field. Examples are astrophysical plasma, solar atmosphere, magnetosphere, 
thermonuclear fusion and liquid metal etc. Their motion described by the MHD equations (see 
[1] for the details), which coupled the Navier-Stokes equations and the Maxwell's equations. 

There are many different approaches to the MHD equations. The wellposedness of these 
equations has been well studied (e.g., cf. [21 [3l HJ [U EJ [7] ) . Numerical methods for compressible 
case were developed in[8l[9l[l0l[IIl[l2l[l3l[Il[l5l[l6l[I71[ll], etc. One can also find similar results 
for the incompressible case in [H [191 l20l \2T\ [22l [23] . The Lagrange invariance of a generalized 
ion magnetic helicity was established for Hall MHD in [23]. Vodak [25] proved the existence 
of a weak, local in time solution to the compressible case on the whole space under special 
assumptions on pressure and entropy. For the historic introduction of the MHD equations, we 
refer the Section 1 of [26] . 

In this paper, we investigate an incompressible case, that is, the MHD equations of incom- 
pressible viscous fluids with finite electrical conductivity (e.g., cf. [27]): 

V-v = 0, (1.1) 

vt + (v • V)v - /io(H X rotH) + -Vp = i/Av, (1.2) 

P 

V-H = 0, (1.3) 

Ut = rot(v X H) + r/AH. (1.4) 

Here v = {u,v,wY is the velocity vector, H = {H^,H^,H^Y is the magnetic field vector, u is 
the viscosity coefficient, /io is the magnetic permeability and r] is the magnetic viscosity. The 
superscript "t" denotes the transpose of a matrix or a vector. 
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Nucci [28 ] found the Lie point symmetries of these equations. Popovych |29[l3m [HT] reduced 
(jl.ip - (jl.4p to ordinary differential equations and to partial differential equations in two and in 
three independent variables by means of Lie symmetries. Then he constructed some classes of 
similarity solutions. 

In order to solve the three-dimensional Navier-Stokes equations, Xu [32] introduced a method 
of imposing asymmetric conditions on the velocity vector with respect to independent variables 
and a method of moving frame which reflect the rotation nature of flow. By these methods, 
he found seven families of non-steady rotating asymmetric solutions with multiple parameter 
functions. In particular, a family of singular solutions may reflect important features of vortices. 
The other solutions are globally analytic with respect to spacial variables. The goal of this work 
is to find similar solutions for the MHD equations (ll.ip - (|1.4p . 

Note that when H = 0, the MHD equations (ll.ip - ()1.4p become the three-dimensional Navier- 
Stokes equations. The factor of magnetic field makes the MHD equations much more complicated 
than the Navier-Stokes equations. For instance, Xu ^2] reduced the Navier-Stokes equations 
under certain conditions to a linear differential equation, and we can only reduce the MHD 
equations to a system of linear differential equations. In particular, we have to use Wronskian 
determinants to solve such systems. Moreover, the magnetic field also naturally leads the ap- 
pearance of Bernoulli equation in our reductions. 

Since the system of MHD equations is a coupled system of Navier-Stokes equations and the 
Maxwell's equations, Xu's methods and exact solutions of the Navier-Stokes equations do give 
us hints of how to solve the MHD equations exactly. We find twelve families of solutions of the 
equations with multiple parameter functions. One of them is a family of singular solutions that 
may refiect important characteristics of vortices. The other solutions are globally analytic with 
respect to the spacial variables, related to polynomials, exponential functions or trigonometric 
functions. Our solutions may help mathematicians to develop the method of proving the well 
known regularity of the equations. They may also help engineers to develop more effective 
algorithms to find physical numeric solutions to practical models. Below we give more technical 
details. 

For convenience, we always assume that all the involved partial derivatives of related func- 
tions always exist and we can change the order of taking partial derivatives. We also use prime 
' to denote the derivative of any one-variable function. 

According to [28], the MHD equations are invariant under the following transformations: 

1. Time translation: 

Ti : t^t + a, (1.5) 

where a is an arbitrary real constant. 

2. Scale transformation: 

Ts : {t, X, V, H,p) ^ (X^t, Ax, Av, A'^H, X~'^p), (1.6) 

where A is any nonzero real constant and x = (x, y, z)*. 

3. Rotations: 

n : {t, X, V, H, p) ^ {t, iix, i?v, RH, p) (1.7) 

for any R G SO{3). 

4. Moving coordinates: 

T4: (t,x,v,H,p) ^ (t,x + F(t),v + F'(t),H,p-px-F"), (1.8) 

where F(t) = {f{t),g{t),h{t)Y is an arbitrary vector valued functions. 

5. Pressure changes: 

Ts: (t,x,v,H,p)^(t,x,v,H,p + 0(t)) (1.9) 



for any function 9 of t. 



The above transformations transform a solution of the MHD equations into other solutions. 
So we only need to find exact solutions modulo these transformations. 
Assuming 



a 



u = ^x + y(j){t,x'' + 



a 



-az + ^lJ{t,x + y ) 



■w 



2 
2 

2 I „,2\ 



H^ = ^y-xi{t,x^+y^) 



-I3z + a{t,x^ + y'^) 



(1.10) 



for some functions a and /3 of t, and two- variable functions (f),ip,^,a, we obtain the following 
solution: 



a 

u = -x 



'^y-xY,am{t){x' + y^r, 



+ yY^a^{t){x^+y^r, 



m=0 
n 



= -x + y^Cm{t){x^ + 



.2\m 



m=0 
n 



m=0 



W 



-az 



+ Y,9mmx'+y'r, 



m=0 



H' = i,y-x^cmmx^ + y'r, (1.11) 



m=0 



,2\rn 



m=0 



and p is given by ()3.86p . where n is any given positive integer, 

k + 2A:2/e-^-* 



a 



(1.12) 



with constants k and / is obtained by solving a related Bernoulli equation, h is an arbitrary real 
constant, and am{t)-,Cm{t)i fmit)-,9m{t) are determined via various Wronskian determinants (cf. 
(j3.76p - (j3.85p ). This is an example of asymmetric approach. 

By Xu's moving-frame approach, we obtain the following solution of the MHD equations 
([nD-dlll): 



u 



, a" HQob , 
'2^ + ^^^-"^ 



6i/ cos a 



4a' 



V =a X — { — - H — - — -)y 

2a' Aa' x cos a + y sin a 

a" fioab 
w =2( — - H ]z 



fj,oab 



X cos a + y sin a 
sin a 



+ Qu sin a 
6z/ cos a 



—X sin a + y cos a 
(xcosa + ysinQ;)^ 
X sin a + y cos a 



{x cos a + y sin a)^ ' 



(1.13) 



4a' 

' = — [{a cos a — b sin a)x + {—a sin a + b cos a)y] sin a, 
= [{a cos a — b sin a)x + {a sin a + b cos a)y] cos a, 



(1.14) 



-6z 



and 



a" /iofl&^/ 
2^ 1^ 
a'' fJ-oab^i 
2^ ^ 1^ 



2, (2; cos a + y sin a) 



4a' 



f-u— w , .a" , /^oo?'s2 / ,.2A-xsina + ycosaf 



4a' 



_ a" , /foo^v , r,/ «" , Atoa^N2N^^ 
^4a'+ 4a' ^ +^4a'+ 4a' ^ ^ 2 

-|-(a" + 2a'(-^ + )){x cos a -|- y sin a){—x sin a + y cos a) 
2a' 4a 



,— xsina + ycosa 
-12i/a' h 



12z^^ 



a; cos a -|- y sin a (xcosa + ysina)^/ ' 



(1.15) 



where a is any function of t, and a, b are arbitrary real constants. The above solution may 
partially reflect important features of vortices. 

The paper is organized as follows. In order to get initial feeling of the MHD equations, we 
solve the equations in Section 2 by assuming that the solutions are linear in the spacial variables. 
Asymmetric approaches are presented in Section 3. In Section 4, we obtain five families solutions 
by the moving-frame method. 

2 Linear solutions 

In this section, we find four families of solutions of the MHD equations ()l.ip - (jl.4p . which are 
linear with respect to the spacial variables x, y and z. 
Denote 

Wi=ut + uu^ + vuy + wu, - fio{H\H^ - H]) - H\hI - H^)) - i^Au, (2.1) 

W2 = vt + uv, + vvy + wv, - MH'^iH'^ - Hi) - H\hI - H])) - u/\v, (2.2) 
W:i = wt + uw^ + vwy + WW, - fiQ{H^{Hl - Hi) - H'^{H^ - H^)) - i^Aw. (2.3) 
To solve ()1.2p . we first deal with the compatibility equations 

Wiy = W2^, Wu = W3^, W2z = W^y (2.4) 

and then find p. The equations (jl.4p can be written as: 



H}- 




-vH^)y + {wH^ 


- uH^), 


- i]AHi 


= 0, 


(2.5) 






- wH^), + {uH^ 


- vH% 


-7]AH2 


= 0, 


(2.6) 




{wH^ 


- uH^)^ + {vH^ 


- wH\ 


-r]AHs 


= 0. 


(2.7) 



Let 




H 




(2.i 



where A = (oij(t))3x3 and B = (6ij(t))3x3 are 3x3 matrices whose entries are functions in t. 
Then by and ([0|), ti{A) = tr{B) = 0. Now the equations ([231) become 



(2.9) 



a'i2 - ^21 + (flii + «22)(ai2 - a2i) + 013032 - 023031 

+ fJ-oiihl + fe22)(fel2 - ^21) + ^13632 - ^23^31) = 0, 
a'l3 - '^31 + («11 + «33)(oi3 - 031) + O23O12 - O21O32 

+ fJ-oiihl + fe33)(^13 - ^31) + ^12^23 - ^21^32) = 0, 
023 - '^32 + (022 + 033) (023 - 032) + O13O21 - O12O31 
+ Aio((b22 + fe33)(^23 - ^32) + ^13^21 - ^12^31) = 0. 

Moreover, the equations (jl.4p are equivalent to the following system of ordinary differential 
equations: 

' b'n = 012621 - 021612 + 013631 - 031613, 

b'i2 = aubi2 - O12611 + O12622 - O22&I2 + O13632 - O32613, 



b' 



13 



'21 



22 



23 



2oii6i3 - 2013611 + 022613 - 013622 + 012623 - O23612, 

-■ -O1162I + O21611 + O2262I - O21622 + 023&31 - O31623, 
■ — O1262I + O21612 + O23632 — O32623, 

-■ O11623 - O2361I + 2O22623 - 2O23622 - O1362I + O21613 



(2.10) 



(2.11) 



and 



631 = -2011631 + 2031611 - 022631 + 031622 - 021632 + 032621, 

^32 = O11632 + 032611 - O12631 + 031 612 - 2022632 + 2O32622, 
, ^33 = -O13631 + O31613 - 023632 + O32623. 

We consider the following special cases: 



(2.12) 



Case 1: A = AK B = BK 



Easily see that the system (j2.9p holds naturally. By (j2.10p - (j2.12p . we find that i? is a constant 
matrix and A is determined by following equation: 



bi2 


-bii + b22 


^23 —^32 


-bi3 


26i3 


b23 


-26ii - 622 613 


-bi2 


^23 


bl3 


— bl2 b23 


-611 - ^ 



:^22/ 



/oii(t)\ 

Ol2(t) 

ai3{t) 

022 (t) 

V023(t)y 



Note that above system has nontrivial solutions. They span a linear space. Since 

Px = -pWi, py = -pW2, Pz = -pWs, 

one gets 



(2.13) 



(2.14) 



p 



((o'li + J2 + ("'12 + Yl + («33 + Yl 4>') 



333 

-Pi{a'i2 + Y 0''^j"'2j)xy + (o'i3 + ^ aija3j)xz + (023 + ^ a2ja3j)yz), (2.15) 
j=i j=i j=i 



modulo the transformation in (11.9 



Proposition 2.1. Let B he any 3x3 constant symmetric matrix with zero trace and let A he 
a symmetric matrix function with zero trace determined by (|2.13p . Then we have the following 
solutions of the MHD equations (|l.ip - p.4p .' 





fx) 






v = A 


y 















and p is given hy (j2.15p . 



(2.16) 



Case 2: A = A\ B = -BK 

As in case 1, one shows that B is constant and 

6?, + 6?, - 262 



On 



12 ^^^13 -23^(,)^ ^^^^ 6f, 26f3 + 6i3 ^^^^^ 



3612613623 

O12 — , Oi3 



n2 



3612613623 
9it) 9{t) 

T ) O23 — T , 

Ol3 O23 



where g is an arbitrary function of t. 



(2.17) 
(2.18) 



Moreover, 

3 3 

P = (Ki + ^a?i + 2/xo(fe?2 + ^?3))^' + K2 + IZ«2, + 2/f^o(&?2 + bL))^' 



3 3 

+ («'i3 + X] °^3j + 2/io(&i3 + ^23))^^) - /o((a'i2 + X] "ii"2i - '2l^obi3b23)xy 

3=1 i=i 
3 3 
+ («'i3 + X] aij«3j + 2fiobi2b23)xz + {a'23 + X a2ja3i - 2fj.Qbi2bi3)yz) 

j=i 3=1 

modulo the transformation in (11.91). 



(2.19) 



Proposition 2.2. Let B 6e a 3 x 3 constant skew- symmetric matrix and let A he a symmetric 
matrix function with zero trace determined by (|2.17|) and (|2.18|) . Then we have the following 
solutions of the MHD equations (jl.ip - (jl.4p ; 



(2.20) 





(A 























and p is given by (|2.19p . 

Case 5.- A = -A\ B = BK 

One shows that A is constant and 

(^11, ^12, fci3, ^22, ^23)* = exp(Mt)c 
by (j2.9p - (j2.12p . where c is any five-dimensional constant vector and 



M 



(2.21) 



Moreover, 



P 



( ° 


2ai2 


2ai3 








-2ai2 





a23 


2ai2 


«13 


-2ai3 


-023 





-ai3 


«12 





-2ai2 








2a23 


\ -^23 


-ai3 


-ai2 


-2a23 





(ai2 + ai3 


)—- 
2 


(ai2 + a 


1/2 

2 \y 

23^"^ 


(«13 



(2.22) 



-ai3a322;y - ai2a23XZ - a2iai32/^) 
modulo the transformation in (11.91). 



(2.23) 



Proposition 2.3. Let A 6e a 3 x 3 constant skew- symmetric matrix and let B be a symmetric 
matrix function with zero trace determined by ()2.2ip . Then we have the following solutions of 
the MHD equations (fTTT]) - (f01) .- 



(2.24) 

\zj \z J 

and p is given by (j2.23p . 

Case J: A = -A\ B = -BK 





fx) 






v = A 


y 








Kzj 







Similarly, we set 



and 



023 -ai3' 
N = I -023 ai2 
ai3 -0-12 



512,^13,^23)* = exp(iVt)c, 



(2.25) 
(2.26) 



where c is a three-dimensional constant vector. Then we get the following proposition: 



Proposition 2.4. Let A be a 3 x 3 constant skew- symmetric matrix and let B be a skew- 
symmetric matrix function determined by (12.260 . Then we have the following solutions of the 
MHD equations (fTTD - ffLll) ; 



(2.27) 



and 



P = -p((-(a?2 + a?3)+2/Uo(ft?2 + b?3))Y + (-(«?2+«i3) + 2/Uo(6?2 + 4))^ 



+(-(a?3 + 023) + 2/io(6i3 + ^23))y + (-^13032 + '2l^obi3b32)xy 

+ (-012023 + 2fIobi2b23)xZ + (-021013 + 2/^,0621613)7/2). 

3 Asymmetric Approaches 



(2.28) 



In this section, we will solve the MHD equations (|l.ip - (|1.4p by imposing asymmetric assumptions 
on V and H, following |32j . 
Firstly, we set 



' u = f + {/3 - il)z)x, 
V = —f3y, and < 

w = 



(3.1) 



where / and g are functions oft, y and z; if and are functions of t and z; (3 and 7 are functions 
of t. 

By dUD-dl^]), 



Wi = ft + fiP - tpz) - Pyfy + Ipfz - t^ifyy + fzz) " fJ-oilVgy " V>9z) 
+ ((/? - Ipzf +13'- Ipzt - Iplpzz + Vljjzzz - tJ-OipiPzz)x, 

W2 = {0^ + P')y - fj-ogyg - nogyij - ^z)x, 



and 



W3 = ipzil + tp) - i^tpzz - iJ-oigzg + {gzil - ^z) - g^zz)x - ^Pzz{i - ^z)x'^)- 

According to (II. 4p and (12. 4p . we find that 7 = 6/2 is a constant and 

Lft = -ip^Pz + tpz^ + rjipzz, 

gt = -7;iyfy - f) + Pi9 + ygy) - i^zg - i^gz + f^z + /zv + vigyy + gzz), 



fyt - Ipzfy - Pyfyy + Ipfyz " l^ifyyy + fzzy) " fJ'Oi-^ygyy + gy^z - fgyz) = 0, 

fzt + fzf^ - fi^zz - Pyfyz + Ipfzz - ^{fyyz + fzzz) " ^0(^^1/2 " 'fdzz " Ij^Qz + g^z 



(3.2) 
(3.3) 

(3.4) 
(3.5) 



= 0, 
(3.6) 



- Wi^zz + i'zi^zz - i'zzt - i^^zzz + '^'(Pzzzz - fJ-oib - VzVzz + ^V'zzz) = 0. (3.7) 

One shows that the system (j3.6p is impHed by 



ft + {fi- ^Pz)f - Pyfy + i^fz - v{fyy + /^^) - fioi^ygy - (pgz + gi^ - v>z)) = 0. (3.8) 
We assume that 

ip = az + l. (3.9) 

Then 

'4: = ahz + hl (3.10) 

a 

by (|3.5p ^. where a is a constant and /i, / are functions of t. Hence the equations (|3.5p - (|3.7p can 
be written as 

( I' b b 

ft + {P- ah)f - Pyfy + {h{az + I) - -)fz - f^oi^ygy - {az + l)gz - ~ "')9) 

- l^ifyy + fzz) = 0, 

b I' 
9t + i^iyfy + f) - {f^- ah)g - Pygy + {ahz + hl )g^ - af - [az + /)/^ - ri{gyy + g^^) = 0. 

^ Z Qj 

(3.11) 

Now we solve above equations. Firstly we suppose that I is a constant and then use the method 
of separation of variables. Set 

f = y{az + l)s, g = y{az + l)q, (3.12) 

where s and q are functions of t. Substituting (13.12P into (j3.1ip . one gets that s is a constant 
and 

q = exp(2 jip- ah)dt){{2a - b)s j exp(-2 j {(3 - ah)dt)dt + c), (3.13) 
where c is a constant. Furthermore, 

2 2 1 

= -p{{(3' -ah' + {[3-ahf)'^ + {-(3' + (3')^ + —{h' + ah^)a''z^ 

2 2 2a 

-fioq'^a^z'^Y ~ l-'-ol^^-^ - o)xyz) (3.14) 

modulo the transformations in (11.81) and in (11.91). 



Theorem 3.1. For any functions h{t), f3{t) and constants a, b, c, and s, we have the following 
solutions of the MHD equations (jl.ip - (jl.4p .' 

u = asyz + (/3 — ah)x, 

< v = -f3y, (3.15) 
^ w = ahz, 

= ayz exp{2 J iP - ah)dt){{2a - b)s J exp(-2 ah)dt)dt + c) + (^ - a)x, 

< ^2 ^ _b_y^ (3.16) 

^H^ = az 
and p is given by ()3.14p . 

Now we give another solutions of (j3.1ip . Motivated from (j3.12p . we assume that 

(f = Cy' + Dz' + Eyz + Fy + Gz + H, 
\ g = ly^ + Jz^ + Kyz + Ly + Mz + N, 



where C, D, E, F, G, H, I, J, K, L, M and N are functions of t. 
Substituting ()3.17|) into (j3.1ip . we get that 



C'-{(3 + ah)C-^io{^ + a)I = 0, 

D' + {P + ah)D + /io(^ + a) J = 0, 
E' = 0, 

/' 

F' - ahF + ihl )E - a^oL + ^iqIK = 0, 

a 

G' + f3G + 2(hl - -)D + -^xqM + 2^iqIJ = 0, 
a 2 

/' h 
H' + {13- ah)H + {hi - -)G - 2u{C + E) + hqIM + no{- -a)N = 

a ^ 



and 



Q 7 

/' - (3/3 - ah)I + (y - a)C = 0, 

J' -{f3- 3ah)J + - Sa)D = 0, 
K' - 2{f5 - ah)K + 2(6 - 2a)E = 0, 

I' 

L' - (2/3 - ah)L + {h - a)F + {hi - -)K - IE = 0, 

a 

b I' 
M' -{13- 2ah)M + (- - 2a)G + 2{hl - -)J - 21D = 0, 
2 a 

h V 
N'-{(3- ah)N + {-- a)H + {hi - -)M - 2r]{I + J) - IG = 0. 
2 a 



Thus we have that 



(2 +a)A'o 



{G' -{p + ah)G), 



C" - 4/3C" - (/?' + ah' - (3/3 - ah){f3 + ah) - (y - a)(- + a)fio)C = 0, 



and 



J 



1 



{^+a)fio 



{D' + {(3 + ah)D), 



D" + 4ahD' + (/?' + a/i' -{(3- 3ah){P + a/i) - (^ - 3a)(^ + a)^o)^ = 0, 

\k' - 2{f3 - ah)K = -{b- 2a)E, 

1 /' 

L = {F' - ahF + {hi - -)E + ^^olK), 

ajjLQ a 

F" - 2(3F' + {ah{2/3 - ah) - ah' + {b - a)afio)F = fF{t), 
M 



-^{G' + (3G + 2{hl - -)D + 2fiolJ), 



N 



G" + 2ahG' + (/3' - /3(/3 - 2a/i) - ^(^ - 2a)Mo)G' = /G(t) 



1 /' 

-V — (//' + (/? - a/i)i? + {hi - -)G - 2v{G + ^) + /io/M), 

(I - a)iiQ a 



H" + (/3' - ah' -{(3- ah)^ -{-- a)^fio)H = fH{t), 



where 
frit) : 

fcit) : 

and 



I' I' I' 

-{{hi - -)E + hqIK)' + 2(/3 - ah){{hl - -)E + ^qIK) - a^o{{hl - -)K - IE), (3.26) 
a a a 

v v v 

-2{{hl - -)D + hqU)' + 2{I3 - 2ah){{hl - -)D + hqU) + bfio{{hl - -)J - ID), (3.27) 
a a a 



1' I' 

{-{hi - -)G + 2u{C + E)- fiolM)' + {P- ah){{hl - -)G - 2u{C + E) + fiolM) 
a a 



+(J - a)fxo{{hl - -)M - 2ri{I + J) - IG). 

Ob 

To write down the solutions of above systems exphcitly, we assume 

/3 — a/i = 



(3.28) 



(3.29) 



where g is a constant. Under this condition, the equations (|3.20P o. (j3.2iP o. (|3.23P 2 and (j3.24P a 
can be written as the following form via certain changing of variables: 



J^'{t)-,kA{t) = fA{t). 



(3.30) 



Here A: is a constant. We will explain this by solving ([320]) 2 in detail. Since the other equations 
are similar, we will only write down the solutions but omit the solving procedure. 
Set 

C = exp(2 / pdt)C. 



Then (j3.20P o becomes 



C" -{q'-{ 



3b 



a){- + a)no)C = 0. 



(3.31) 
(3.32) 



We write Ac = q"^ 



(y — o)(| + a)/io and define 



ex.p{y/X^ t + 2f pdt) if Ac > 0, 
exp(2 / pdt) if Ac = 0, 

cos(iv^ t) exp(2 / I3dt) if Ac < 0, 



Then 



and 



exp(-VA^t + 2//3dt) if Ac > 0, 

= { iexp(2 / Pdt) if Ac = 0, 

sm{i^/X^ t) exp(2 / /3dt) if Ac < 0. 

C = IciCci + IC2CC2 



E/c,(c^^,-(2/?-g)ec,) 



+ a)fio 



(3.33) 

(3.34) 
(3.35) 

(3.36) 



where Zci and /c2 are real constants. 



Similarly, we set Ad = 5^ + (| — 3a)(| + a)/.fo and define 

exp{^/)^t-2J{|3 + q)dt) ifAD>0, 

Cdi = { exp(-2 /(/? + q)dt) if Ad = 0, 

cos{i^/XE t) ex.p{-2 J {(3 + q)dt) if Ad < 0, 

exp(-VA^ t-2j{l3 + q)dt) if Ad > 0, 

^D2= { t exp(-2 J{P + q)dt) if Ad = 0, 

sm{i^/XB t)exp{-2 f{l3 + q)dt) if Ad < 0. 



(3.37) 



(3.38) 



Then 



and 



J 



E/D,(e;5, + (2/3-^)CD,) 



(I + a)/io 



where //ji and lo^ are real constants. 

Given functions r{t) and we denote their Wronskian determinant by 



W[r{t), s(t)] = det f^^^^ = r{t)s'{t) - r'{t)s{t). 



For notational convenience, we write 



WMt),s{t)]= det (^^ '^^^j^ 



-s{t), 



Set \f = 



Then 



and 



P72[r(t),s(t)]=det( f) =r(t). 



r'(t) 1 



a)afiQ and define 



exp(v^ t + / /?dt) ifAF>0, 
= <( exp(//3dt) if = 0, 

cos(i-v/Xp t) exp(J /3dt) if Ai? < 0, 

' exp(- VAf t + / (5dt) if Af > 0, 
^F2 = s texp(j' (5dt) if Ai? = 0, 

sin(iVXp t)exp(/ (3dt) if Xp < 0. 

' ' Jo W[^fAs),^F2{s)\ 



L 



Jo 



'W,[Cf,{s),Cf2{s)] 



W[^FAt),CF2{t)] 
fF{s)ds))), 



W[^F,is),^F2{s)] 

where Ifi and If2 are real constants. In this expression, we have use the fact: 



2q 



where k and E are real constant 

b(b 
2\2 



Let Ac = 9^ + I (I — 2a) /io and denote 



exp(VA^t-/(/3-g)dt) 
Cg, = { exp(-/(/3-g)dt) 



if Ag > 0, 
if Ag = 0, 



^G2 



cos{^^/XG t) exp(- - q)dt) if Ag < 0, 

exp(- VA^ t- J{(3- q)dt) if Ag > 0, 

t exp(- J {(3- q)dt) if Ag = 0, 

sin(iVA^ t)exp(- /(/? - (?)dt) if Ag < 0. 



Then 



and 



2 / (R-q-l'Vj^. I 



' ' ' Jo W[^Gi{s),^G2{s)\ 

+^^^w^[eG.(t),eG.(t)]^^^^^^> ^'-''^ 



where /gi and I'sal constants. 

^2 



Set Xh = + (| — a)/io and define 



ex.p{^/XH t) if Ah > 0, 

Ch^ = { 1 if A// = 0, (3.53) 

cos{i\/XH t) if Xh < 0, 

exp(— a/Xh" t) if A// > 0, 

Ch2 = { t if A// = 0, (3.54) 

sm{^^/XH t) if Xh < 0. 



Then 



and 



Jo W[iHAs),iH2[s)\ 



(3.55) 



where -E, ///j and Ih2 are real constants, C, G and M are given by ()3.35|) . (|3.51|) and (j3.52p . 
Furthermore, by (j2.14p . one has 

p = -pf^x^ + {H' + qH+ ^^~'^^^~^' G + nolM - 2u{C + D))x 



-^J^z^ - fioJMz^ + {13' + {[3- qf - fio{2JN + M^))z^ 

+ (-//oMiV + -(/?'/ + {(3- q)l' - I") + ^^{{(3 - q)l - l'))z 
a a 

-(^ - a)noKxyz - p^Ixy^ "-^t^^Lxy - (^ - a)^iQJxz^ «)/^oAfxz 

-fioKJyz^ - fioKIy^z - y (2/J + K'^)y'^z'^ - pq{JL + KM)yz'^ 



-po{KL + IM)y' - noiKN + LM)yz ) (3.57) 
modulo the transformation in (jl.Op . 



Theorem 3.2. For any functions f3{t), l{t) and real constants a, b, q, E, we have the following 
solutions of the MHD equations (jl.ip - (jl.4p ; 



' u = Cy^ + Dz^ + Eyz + Fy + Gz + H + qx, 
V = -Py, 

{(3-q)l-l' 



w 



= {f3- q)z + 



(3.58) 



= ly'^ + + Kyz + Ly + Mz + N + a)x, 



H'^ = az + l 



(3.59) 



and p is given by (j3.57p . Here the functions C, D, F, G, H , I, J, K, L, M and N are given by 
(I335D. (I539D. (i06D. (133X1). (I335I). (I336D. (13301), (i08D. (iOTD. (1332!) and (13361). respectively. 



We denote 
and assume 



2 , 2 

= x + y 



= -j^x + y(l){t,w), 



a 



v = ^y- x(l){t,w), 
^ w = —az + tp(t, vj), 
Recall that (jl.4p is just 

( Hi = iuH' - vH% 
Hj = {vH^ - wH'^)z 



-x + yi{t,w) 

-y - xiit.w) 
H^ = -I3z + a{t,w) 



H' 



H' 



{wH^ -uH^), + riAH\ 
{uH^ -vH^)^ + r]AH^, 



, Hf = {wH^ - uH^)^ - {vH^ - wH\ + T]AH^. 

Substituting (|3.6ip into (j3.62p ^. one gets that /? = 6 is a constant and 

= -aro^ro + bzucj)^ + 4r/(tz7^)„^. 

One shows that (|3.62P n is equivalent to (j3.63p . and (|3.62P n implies that 

at = b{zuip)^ - avaa^ + Ari{a^ + wa^-ai)- 

Moreover, by (j2.4p . we have 

J (ro(^)to + a{w{wcl))^)^ - Av{w{zD(l))^^)^ + b^Q{zD{Tui)^)^ = 0, 
I il)t- ail) + awip^ - Au{'ip^ + zutp^^) + bfiQwa^ - bfiQU = 0. 



(3.60) 



(3.61) 



(3.62) 



(3.63) 



(3.64) 



(3.65) 



Given any n E N, set 

n n n n 

<P=Y1 ^ = E cm(i)^", ^=Y1 /-(*)^™' ^=Y1 9mit)w"'- (3.66) 

m=0 m=0 m=0 m=0 

Substituting (j3.66p into (|3.63p - ()3.65p and considering the coefficients of w"^, we obtain that 



am,t + [m + l)aam + bno{m + l)cm = 4z^(m + l)(m + 2)0^+1, 
Cm,t - hmam + amcm = ^ri{m + l)(m + 2)cm+i 



(3.67) 



and 

J fm,t = -amfm. + b{m + l)gm + 4:r]{m + iffm+i, 
1 9m,t = -{m - l)agm - b{m - l)^o/m + 4z^("i + lYgm+i, 



2 (3.68) 



where a_k = Un +k = C-k = Cn+k = f~k = fn+k = 9~k = 9n+k = for /c > 0. These systems are 
similar to ()3.20p . If we set 

k + 2/c2/e"^'* 

a = — — ^-TT (3.69) 

for any constants k and then we can transform them into certain second order ordinary 
differential equations with constant coefficients, as we did in (13.20p . 
Observe that 

= --, TTTT («m,t + (w- + 1)00.^ - ^vivTL + l)(?n + 2)a.m+i). (3.70) 

Substituting p.70p into (I3.67p . one has 

+ (2m + l)aa^ + (m + X){rno? + q' + }? [iQm)am. = Fm{t), (3.71) 

where 

Fm = 4:{m +l){m + 2){miyaam+i + i^a^.+i + ?/c,m+i). (3.72) 

Set 

am = am exp(-^ J {2m + l)Qdt) = (e(™+^)^* - 2A;/e(™-5)'=*)a™. (3.73) 
Then (j3.7ip becomes 

2iji 4-1 1 

+ ((m + l)(ma2 + a' + 6Vo"i) — a' - -{2m + lfa^)dm = Fm{t). (3.74) 

By 

2m +11 1 

{m + l){ma^ + a' + b'^fiom) — a' - -{2m + ifa^ = m{m + l)b'^fio - (3.75) 

is a constant. Let Am = —m{m + 1)6^ /-to + \k? and define 

exp(VA;;; t)(e(™+5)'^* - 2A:/e("^-5)*^*) if > 0, 

Cm,i = -i e(™+^)'=* - 2A:Ze('"-^)^* if A„ = 0, (3.76) 
cos(iv^ t)(e('"+5)'=* - 2A:Ze(™-5)'=*) if A^ < 0, 

' exp(- Va;^^ t)(e('"+i)'=* - 2A:/e(™-^)'=*) if A^ > 0, 

Cm,2 = < t(e("'+5)'^* - 2/k/e("^-5)*^*) if Am = 0, (3.77) 
_ sin(iVA;;; t)(e(™+^)'=* - 2A;/e(™-3)'=*) if A^ < 0. 
Then we can get am and Cm by the following recursive formulae 

V^A _L ^i[Cm,l(s),Cm,2(g)] ^ , . 



hi{){m 



1 , r ^,[Cm,l(^),C»x,2(g)] ^ 

„ + 1) IWi^™-.- + iy[Cm,l(.), UM\ 



+Cmj'( 



W"j[Cm,l(0,Cm,2(^)] 
W^[Cm,l(t),Cm,2(t)] 



Fm(t) + {m^\)a{l 



(s)ds) 

'o W^[Cm,l(s),Cm,2(s)] 



+4i/(m + l)(m + 2)Cm+i,i(Wi,i + / 

JO 



W;?-[Cm+l,l(g),Cm+l,2(g)] 
W^[Cm+l,l(s),Cm+l,2(s)] 



Fm+i(s)ds)), (3.79) 



where Sind 1^,2 are real constants. 
Given 

Gm{t) = 4(m + lf{{m - l)ar]fm+i + vfL+i + ^9ni+i] 



(3.80) 



the solutions of system (|3.68p can be written down similarly as those of (|3.78p and (j3.79p . 
Set 



and denote 



eMV^ t)(e('"+i)'=* - 2kle^^-\^^') if r„ > 0, 

= { e('"+5)'=* - 2A;Ze('"-^)'=* if r„ = 0, 

cos(iV^ t)(e("^+3)*^* - 2fcZe("*-^)'=*) if r„ < 0, 

exp(-^ f)(e('"+5)'=* - 2A;/e("^-5)'^*) if > 0, 

?m,2 = <; t(e("^+5)'=* - 2fcZe("-5)'=*) if r„ = 0, 

sin(f^ t)(e(™+5)'=* - 2A:/e('"-5)*^*) if r„ < 0. 



Then 



and 



Jm — / ^^m,j\qrn,j + / TTTT 7~TT'^m\S QS) 



1 ^ f 



•s),?m,2(s)] 
* W"j[?m,l(s),W,2(s)] 



(3.81) 



(3.82) 



(3.83) 



(3.84) 



W^[<?m,l(s),W,2(s)] 



W^j[';'m,l(^),<?m,2(0] 
VF[w,l(t),?m,2(0] 



Gm{t) + ma{qm,j + / 

JO 



Gm{s)ds) 

W^[?m,l(s),?m,2(s)] 



Gm(s)ds)) 



Jo 



VFj[w+i,i(s), Wi,2(s)] 



l^[Wi,i(s) 5 W+1,2 
where qmi and 5^,2 are real constants. The function p can be written as 



Gra+iis)ds) 



(3.85) 



P = -pi 



a' + 



+ -(2a^ — a')2:^ + bfiQza ~ 2 J 4'^'^'^ ~ + j vj^S^i^dw) 



(3.86) 



modulo the transformation in (jl.Op . 

Theorem 3.3. Given any constant h, positive integer n and the function 

k + 2kHe~^^ 



a 



-1 + 2/cZe-^* 

with constants k and I, we have the following solutions of the MHD equations (jl.ip - (jl.4p .' 



(3.87) 



Q 

U = —X 

2 



+ y^am{t){x^ + y'^Y 



a 
T 



m=0 
n 



V = -y - am{t){x'^ + y 

m=0 



2\m 



H^ = -x + yY^ Cm{t)ix'^ + 



m=0 
n 



w = —az 



+ Y.g^{t){x^ + y^r, 



m=0 



H' = -y-xY, cmit)ix' + y^r , (3.88) 

m=0 
n 

H^ = _bz+^fUt){x^+y'r. 



m=0 



and p is given by p.86p . Here the functions a^, /,„ and gm are given by (|3.78p . ()3.79p . 
(|3.84p and (|3.85p . respectively. 



4 Moving-Frame Approach 



In this section, we use the moving-frame method to find six famihes exact solutions of the MHD 
equations (frT]l - (f01) . 

Let a and /3 be two functions in t. We denote 

(cos a sin a cos /3 sin a sin /3\ 
— sino cosacos/3 cosasin/3 . (4-1) 
— sin /? cos /? / 

Then 

(cos OL — sin a 

sin a cos /3 cos a cos /3 — sin/3| . (4.2) 
sin a sin /? cos a sin /? cos /? 

Following section 3 of [32], we define the moving-frames: 

^=in\=T[ii'\, v=fv)=rft;V (J|=r|y|. (4.3) 

\n^j \Hy \wj \wj 

Write 

\/ = [dx ,dy , dz) and to* M = V x M. (4.4) 

Then we have 

Lemma 4.1. We have the following equations: 

HxT^t H=T{Hxrot H), (4.5) 

fot(v xH) = T{rot {vx H)). (4.6) 
Proof. It follows from straightforward calculation. □ 



We write 



Ri = Ut + a'{yUx-XUy-V) + {ZUx-XUz-W)f3's\xia 

+ {ZUy - yUz)fi' COS a + UUx + VUy + WUz - uA{U), (4.7) 

i?2 = Vt + a'{yVx-XVy+U) + iZVx-XVz)P'sma 

+ {ZVy - yVz - W)P'cosa+UVx + VVy + WVz - i^A(V), (4.8) 

i?3 = )^^ + a'{yWx-XWy) + {ZWx-XWz+ll)P'sma 

+{ZWy - yWz + V)P' cos a +UWx + VWy + WWz - z^A(W), (4.9) 

GT-ZS /'7_/2 \ 'l_/3/'7_/l '1-/2 \ 

1 = n {ttx - rty) - n [nz - nx)-, 

G2 = n^{nl,-nl)-nHnl-nl>), (4.10) 

3 = At [ttz - rtx) - n [ny - nz), 
= Rj - fioGj j = 1, 2, 3, (4.11) 

= nl + a'{yn\; - xn]) - n^) + {zh\ - xnl - n^)p' sma 
+{znl; - yH\)i3' cos a - {un^ - vn% + {wn^ - un^) z - v^'h\ (4.12) 



+{zn^ - yn% + n^)p' cos a - {wn^ - un^)x + {vn^ - wn^)y - tjAh^ 

Then the MHD equations (fTTIl - lfLi]) become 

1 1 1 

^1 + -Px = 0> ^2 + -Py = 0, ^>3 + -p^ = 0, 
P P P 

= 0, j = 1, 2, 3, 

f Ux +Vy + Wz =0, 

Instead of solving (j4.15p . we wih first deal with the compatibility conditions: 
Then we find p via (j4.15p . 



Now we suppose that 




ay — Z(3' sin a, 






= g-2iX 




and < 




= <y + iy. 











where (p, ip, a and h are functions in X and t, and g, 7 and ^ are functions in t. 
Substituting KW\ into (HT6]l - (liT8| ). we get that 



e' = o, 



-275 



and 



(pt - 2jX(l)x +l4>- '^4>xx - ipP' cos a + HQ{^,a + ga^ - 2Xiax) 

= ((/3')^sinacosQ + a'7-a")^ + 6'i, 
il)y - 2-fXipx - vipxx + 0/3' cos Q + noiih + ghx - 2Xihx) 

= —{{13' sin a)' + a' j3' cos a + a' (3' cos a — 7'/?' sinQ)A' + 62, 
at - 2jXax - ja - rjaxx - h(3' cos a + {^4> - gcpx + 2X^(j)x) 
= 3a' X^ — a'g, 

Tit - 2'yXhx - 7/1 - ^hxx + crP' cos a + {^,1^ - gipx + 2X£^ipx) 
= SXP'^ sin a — gj3' sin a, 

where and 62 are any functions of t. 
Set 

(/5 = y /3'cosadt and /i = exp(4 j 7dt). 
We have g = cj ^fji by (j4.20|) . where c is a constant. 



Change variables by 



/A 


/ 






a 




\hj 


V 



_ 1 _ 1 

fi 4 COS If i-i * sin If 

—fi 4 sin If i cos V9 

i_ 1 
^4 cos (p fi-i sin if 

—fiisimp fiicosif/ 



\ 

a 
\hj 



Then ()4.2ip becomes 



u - - c 

o-j - —Xox - Wxx + (px + 2— = + /32, 

2/i ^ fi ^ 

fj.' " " ^-c- 2^- 
ht - —Xhx - rjhxx + '^l^x + -^Xipx = q^X + /34, 

2// ^ H y/JI 



where 



qi = sin a cos a + a'7 — a") cos ip — ((/?' sin a)' + {a' — j')f3' cos a) sin tp), 

q2 = (~((/?')^ sin Q cos a + Q'7 — a") sin — ((/?' sin a)' + {a' — 7')/?' cos q) cos tp), 
qs = fj.~^ (3a'^ cos p + 3/3'^ sin a sin , 
, f/4 = /x~ 4 (— 3q;'^ sin p + 3/3'^ sin q cos p) 



and 



Pi = /-f 4 (^1 cos p — 62 sin 93), 
P2 = f^i-^i-Si sin 93 - 6^2 cos 93), 
Ps = fi~i{—a'gcosp — P'^ sin a sine/?), 
^ Pi = fi~^ {a g sin p — sin a cos 99). 
Let n be a positive integer. We write 



m=0 m=0 m,=0 m=0 

Substituting above expressions into ()4.24p . we obtain that 

am,t - {m + l)—am - (,fioy/Jl{2m + 1)6^ = Am{t), 
A* 

bm,t - {m+ l) — bm + (2m + 3)— a^, = 
2/i ^ 



where 



and 



Am = -cfj,o{m + l)bm.+i + (m + l)(m + 2)1/0^+2 + + 5m,oPi, 

Bm = -{m + l)am,+i + {m + l)(m + 2)r]bm+2 + 5m,ig2 + '5m,o/32 



Cm,t - (m + 1) — Cm - ^Ho^/Jl{2m + = Cm(t), 

c^m,t - (m + I)— dm + (2m + 3)— = -Dm(t), 
2// ^ 



where 



-CHo{m + l)drn+i + (m + l)(m + 2)vCrn+2 + Sm,iq3 + ^mfiPs, 
-(m + l)Cm+l + {ni+ l)(m + 2)l]dm+2 + Sm,iq4 + 5m,o/34- 



As we did before, let 



H = I 



exp( y"(2/c + — 



,)dt) 



(4.31) 



(4.32) 



+ 2/26-^=*' 

for any constants k, li and l2- Under this condition, the equations (j4.28p and (j4.30p may be 
solved as second order constant-coefficient ordinary differential equations via certain changing 
of variables. We only give the results but omit the procedures. 



Write Xm = jk^ - (2m + l)(2m + 3)^>o and define 



Cm.l 



■>m,2 



jj. exp{^/X^, t) if Am > 0, 

jl if Am, = 0, 

IICOs{iy/Xm t) if Am < 0, 

fi exp(- VA^ t) if Am > 0, 

tfl if Am = 0, 

fism{iy/Xm t) if Am < 0. 



In addition, we denote 

Pn 
Qrr 



A' — A 

C„, — Cn 



(m + l)Am + '^Ato(2m + l)^Bm, 



[m + l)Cm + S.fJ'oi^m + l)^Dm. 



Then we get the following recursive formulae: 

2 



(s),?m,2(s)] 



Pm(s)ds), 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



\ Vf,' .+ r ^,[w(^).^m,2(^)] 



W[c,m,l{t),(,m,2{t)] 2/i Jo («) , ?m,2 (s)] 



(4.37) 



Cm= y<im,j{qmj+ -T^, ^QmS)ds, 4.38 

^ Jo VF?m,l(s),S-m,2(s 



(5m(s)ds) 



and 

, _ 1 Sr i < ( . ^j[?m,l(g),^m,2(g)] 

(2m + l)^^oV^^V^'"'^'^^"^'''^yo W^[w(^),?r.,2(.)] 
+ W(#4|^QUt)-(m + l)|:(,m.+ r#^^#^Q™(.)d.)) 

W/[?m,l(t),?m,2(i)J 2// Jo l^[?m,l(s), W,2(s)J 

(4.39) 

where gm,i and gm,2 are real constants. 



Theorem 4.2. Let a and (3 he arbitrary functions of t, and let c and ^ he arbitrary real con- 
stants. The functions ip, n, cj), a, ■0 and h are given in (j4.22p . (j4.32p . ()4.27|) and (j4.36p - (j4.39|) . 
respectively. Then we have the following solutions of the MHD equations (ll.ip - ()1.4p ; 

u = — X-^ cos a — y{a' cos a + — sin a) — Z(3' sin a cos a — ^u"* (0cos — ip sin (p) sin a, 
/i 4/i 

V = — X-^ sin Q cos (3 + 3^(— a' sin a + ^ cos a) cos B + Zi—d' sin cos /? + -7— sin /?) 

/i ^' V 4/x 4// ^ 

+ /i~ 4 ((cos a cos /? cos 99 + sin /3 sin ip)^ + (cos a cos /? sin — sin /3 cos p)'^^^ 

w = — X sin Q sin 8 + 3^(— a' sin a + — cos a) sin /3 + Zi—d' sin sin /? + — cos S) 

[I \ 4fi ^ ^ 4/i ^ 



+ /U 4 ((cos a sin /3 cos 99 — cos Psm(p)(f> + (cos a sin /? sin 99 + cos /3 cos ■;/'), 



= cos a — 2X£^ cos a — sin a — sinQ((7 cos 93 + hsimp), 
v 

-ff = sin a cos /3 — 2^%"^ sin a cos /3 + 3^^ cos a cos /? — ZS^ sin /3 

+ fi^i ((cos a cos (3 cos + sin (3 sin y?)*? + (cos a cos /? sin (yj — sin/? cos ip)h), 

= sin Q sin /3 — 2XS^ sin a sin (3 + cos a cos /? + cos f3 
^//^ 

+ ((cos a sin (3 cos (/? — cos (3 sin (^)(T + (cos a sin /? sin p + cos /3 cos (p)h), 



(4.40) 



(4.41) 



- sin a)2)— - {a" + ^a' - f3' sin ap' cos a)A'3^ 
2 2/x 



-((/3'sinQ)' + Q'/3'cosa + — /3'sinQ)Af^-3^Za'/3'sina-2Q' /" ^d^- 

-2/3' sin a ^ V'd;if - Aio(^(<T2 + /i^) + ^^^3^ + (4.42) 
where (p, ip, a and h are given by ()4.23p . 



Set 

^ U = anX + a^y + ai^Z - QuX~^, 
V = a2iX + a22y + a2-iZ - Gi^yX'^, 
W = 031^ + 0323^ + 033^:, 

Substituting (|4.43p into (j4.18p . one gets that 



'n^ = bnX + bi2y + bnZ, 

re = b2ix + b22y + 623^, 
, n' = bsix + bs2y + bssz. 



an = a22 = 7) 021 = -ai2 = a', (331 = -ai3 = ^'sina, 023 



032 



and 



2a" + 4a'7 = ^0(^31(^32 - ^23) + ^33 (^21 - ^12) + b32(bi3 - ^31)), 
^22(^13 - ^31) = -^23(^21 - bi2) - 621 (&32 - ^23), 
^11(^32 - ^23) = ^12(^32 - ^23) - ^13(^21 - ^12)- 
By ()4.16p . we have that 



b'u^ + b'i2y + h'l^Z + bii{a'y + Z(3' sin a) + 612(73^ - 61/3^^-^) 
- 2(A'/?'sina + 3^/5'cosa + 7Z)6i3 - + &uX~'^){biiX + bi2y + bi^Z) 
+ bii{-iX -a'y - Z(3'sina-QvX-^) = Q, 



(4.43) 



-/3'cosa (4.44) 



(4.45) 



(4.46) 



h'^^X + b'^^y + ^23^ + hiilX - 61/ Af-^) + h22{-a'X + Z 13' cos a) 
- 2{X[3' sin a + yi3' cos a + 7^)623 + (7 + 6z.;f-2)(6i2'^' + ^223^ + 623^) 
+ h22{a'X + 73^ - Z()' cos a - QvyX'^) = (4.47) 



and 



b'^^X + 6^523^ + 6'33^ + 63i(7;f - 6z.^-^) + 632(73^ - QvyX-'^) 
- b^^{XI3' sino + y0 cos a) + 2{b2iX + 6223^)/?' cos a 

+ 2-f{b-iiX + b^2y + h-iZ)-bz^{X 13' sina + y 13' cos a + 2-fZ) = Q. (4.48) 
One shows that H 7^ if and only if /? is a constant. In this case, 

hi = bi2 = bi3 = 623 = hi = ^32 = (4.49) 

and 

6'22 = &21 = 0. (4.50) 

Hence, one gets 

7 = -—-^^^. (4.51) 

' 2a' 4a' ^ ' 

Theorem 4.3. For any function a{t) and constants a, b, we have a solution of the MHD 
equations (fLT]) - (fLl|) .- 

a" fJ-oab . 61/ cos a —xs'ma + ycosa 

u = - {— + -—)x-ay ■ : + Qvsma- ■ 

la 4a xcosa + ysvua (xcosa + ysma)^ 

, a" liQob 6z^sina sin a + 7/ cos a /, 

V =a X — { — -H — - — -)y dz^cosot l^.ozj 

2a' Aa' xcosa + ysma (rr cos a + y sm aj^ 

a" noab 

' = — [{a cos a — b sin a)x + (—a sin a + b cos a)y] sin a, 

= [{a cos a — b sin a)x + {a sin a + b cos a)y] cos a, (4.53) 
= -bz 



and 



f,,a" Uoob . a" UQab r, . ^ (xcosa + ysina 

a" ^J'Gah , a" fiQab 2 ,,2s (— 2;sina + y cosa)^ 

+ + W + ) 2 

_ a" , [J;oab , a" Uoo^ 2n^^ 
^4a' + 4a' ^ + 4a' + 4a' ^ ^ 2 

+(a" + 2a'( — 7 + ))(a: cos a + y sin a)(— x sin a + y cos a) 
2a! 4a 

, — X sin a + y cos a 12z^^ \ .x 

-12i/a' ^— + - : -.]. 4.54 

xcosa + ysma (xcosa + ysmaj^/ 



Let ai, /?! and 7 be functions of t and a, 6 be real constants. Set 

r e°i^+^i^ - ae"("i^+^i^) if r = 0, 
- \ sin(ai3^ + /3iZ) if r = 1; ^^'^^^ 

r e°i^+^i^ + ae-("i^+^i^) if r = 0, . , 

1 cos(ai3^ + /3iZ) ifr = l; ^^'^^^ 



e-fX-be-"/^ ifs = 0, J , r e^^ + be-^^ if s = 0, 
sm(7;^) if . = 1 V. = I ^^^(^ .f ^ ^ ^_ (4.57) 

We assume that 

-{aia + /3iT)Cr0, 
+ I <T-f^ri^s I , (4.58) 





























-{ 











/-(aiCJi + /3lTi)Cr0A 

+ fTnCrV's , (4.59) 

where A = (aij)3x3 and B = (6ij)3x3 are 3x3 matrices whose entries are functions of t; the 
functions a, t, ai and ri are also functions of t; the numbers r and s may be or 1. 
Now 

Ri = {a'li + afi + ai2(a2i — a) + 013(031 — /3'sina) — 021a' — a^i/S' sina)X 

+(a'i2 + aii(ai2 + a') + 012022 + oi3(«32 - /3'cosa) - 0220' - 032/?' sin a) 3^ 
+(a'i3 + aii(ai3 + /J' sin a) + 012(023 + /3'cosq;) + 013033 - 023a' - a^sP' sin a)Z 

-i-lYiaia + PiT)[a[y + P[Z + ai((o2i - a')X + 0223^ + (023 + /?' cos a)Z) 
+/3i((a3i - /3'sina)^ + (032 - /3' cos a)y + a33Z)](,r(ps 
+7[(«i2 - Oi')a + (oi3 - /3'sina)r]^rV's 

-[(aid + /3it)' + on (aid + ^r) + z.(aia + (3^r)i{-iyj'' + (-!)''(«? + /3?))]Cr</'. 
-{aia + /3ir)[7'A' + 7(011;!^ + (012 + a')y + (013 + P' sma)Z)]Cri^s, (4.60) 

-R2 = (021 + 021011 + 022(021 - a') + «23(«3i - /3' sin a) + oiia' - 031/3' cos a)A:' 
+(022 + 021(012 + a') + O22 + 023(032 - P'cosa) +ai2a' -032/?' cosa)>' 
+(023 + 021(013 + /3'sina) + 022(023 + (3' cos a) + 023033 + 013a' - 033/?' cos a) 2^ 

+ {6r,l + 405^,0)0-72 (aicr + (3iT)CrCr 

+(_l)V7[7';f + liauX + (012 + a')y + (013 + P' sin a)Z)]er.</>s 

+ [(cJ7)' + 022CT7 + (023 - /?' cos a)r7 - i/(J7((-l)"72 + (-l)''(a? + /3?))]erV's 

-(021 + a){aia + /3ir)Cr</'s 

+cj7[(a'i3^ + P[Z) + ai((o2i - a')X + 0223^ + (023 + P' cos a)Z) 

+/3i((o3i - P' sin a);f + (032 - P' cos a)y + 033^)]CrV'., (4.61) 

-^3 = (03;^ + 031011 + 032(021 - a') + 033(031 - /3' sin a) + On + 021/3' cos a)A' 
+(032 + 031(012 + a') + 032022 + 033(032 - /3'cosa) + 012 + 022/3' cos 0)3^ 
+(033 +031(013 + /3'sina) +032(023 + /3'cosa) + + 013/?' sin a + 023/3' cos a)Z 
+((5r,i + 4o(5r,o)T7^(aio- + PlT)irQr 

+ (-l)V7[7';f + 7(oii;f + (012 + a')y + (013 + P' sina)Z)]Cr<^s 

+ [(t7)' + (032 + /3'cosa)CT7 + 033r7 - uT-i{{-iy-f'^ + {-lY{a\ + /3i2))]C^'0. 

-(031 +/5'sina)(ai(j + /3ir)Cr<^s 

+Tl[{a'^y + /3i^) + ai ((021 - a')X + 0223^ + (023 + P' cos a)Z) 

+/?i((a3i - /?' sin a)A' + (032 - P' cos a)y + 033^)]CrV's, (4.62) 

Gi = (621-6l2)(62lA' + 6223^ + 623^) + (&31-fel3)(fe31^ + fe323^ + ^33^) 

+7[(-l)^7'(^? + rl) + (-l)^(aiai + Pinf]erCt>s^s 

+ [{b2iX + 6223^ + 623^)((-l)Vi72 + (-l)^ai(aiai + /Jin)) 

Hb^iX + 6323^ + 633^)((-l)Vi72 + (-l)^/3i(aiai + PlTi))]ir(l>s 

+7ki(^2i - 612) + Ti(63i - bi^)%i^s, (4.63) 



G2 = {h2-h2^){h-iiX + h2y + h^zZ)-{h2i-h2){hiiX + hi2y + hi^Z) 
+{axai + /3iri)((-l)Vi72 + (-l)^ai(aicji + f3iTi))irCAl 
+7^ri(aiTi - f3iai)irC,rifl 

-{huX + hi2y + 6i3^)((-l)Vi72 + {-lyaiiaiai + l3iTi))ir'Ps 

+ {hz2 - &23)n7CrV's - (&21 " ^12)(aiO-l + /3iri)Cr(/>s 

+(631^ + + b33Z){aiTi - Piai)jCri^s, (4.64) 

Ga = {bl3-b3i){biiX + bi2y + bi3Z)-{b32-b23){b2lX + b22y + b23Z) 

+{aiai + /3iri)((-l)Vi72 + {-ly pl{aiai + PiTi))Ur<Pl 
-{bnX + buy + 6i3^)((-l)Vi72 + {-lypiia.ai + f3in))ir<Ps 

-{b32 - 623)0-l7CrV's " (&13 " &3l)(aiO-l + /3iri)Cr(/>s 

-(621^ + b22y + b23Z){aiTi - /3i(Ji)7Cr.V's. (4.65) 

We assume that aicJi + PiTi / and aia + f3iT / 0. Then the coefficients of Cr4'r dy^i and 
9x^2 yield 

diTi = (4.66) 
For simpUcity, we consider the special case in which 

021 = — «12 = ct' , 031 = — ai3 = /3'sina, 023 = 032 and B = . (4.67) 

Furthermore, the nonlinear terms in ^ly, $12, ^2X^ ^2Z, ^3X and <I>33; hint us that it is 
convenient to suppose 

^2 = (_l)r+.+l(„2^^2)_ (4_gg) 

Thus the coefficients of Sari's in ^2y and ^2Z suggest that 

aiT = (3ia. (4.69) 

Hence we can write 

ai = ga, n = gr. (4.70) 

As Xu did in [32], we get that 

1"^ = ^"; and 1" = ^"; (4.71) 

for some functions if and fi. Moreover, 

— Oil (a')2 — a'a" + /?' sin a (/?' sin a)' + 2a'/3' sin a /?' cos a 
= {a'y + iP'sinay ' ^^■^^> 

—Oil (/?' sin a)^ + a'a" — /3'sina (/?' sin a)' — 2a'/?' sin a cos a 
= (a')2 + (/3'sina)2 

and 

— oiia'/?' sin a — a"/?' sin a — a'/3" sin a — 2(a')^/3' cos a + (/3')^ sin^ a cos a , , 
= (a')2 + (/3'sina)2 • ^^■^^> 

The coefficients of Cri^s and y in ^'i = show that 

bii = 612 = bi3 = 0. (4.75) 



Substituting i^M^-i^T^ into (glSD, one gets that 

y (ff(aia + /3iT)y 
«ii = = —7 , ^ N . (4-76) 

^22(023 + 032) = ^33(022 - 033), (4.77) 

{boo = 2623/3' cos a, , , 

f , (4.78) 

623 = -2622/? cos a, 

^22 = -^(aia'i - + a\a22 - PiPfa^s - 2qi/3i/3' cos a), 

^23 = ^((ai^i)' + ai/?i(a22 + 033) + (^23 + /?' cos a)r^75r + (023 - /?' cos 0)0-^79) 
and 

-a{gajy + r^grj)' + a22'7^7<7 — a^sr'^jg + 2aTjg(3' cos a 



^^"' + "') ' (4.80) 

-T{gajy + a{gTjy - auarjg + (023 + /3' cos a)T'^jg + (023 - cos 0)0-^75 



^23 



7(cr2 ^ ^2^) 

One shows that the above systems (j4.77p - (j4.80p are compatible with each other if and only if 
a"P' sin a - a' (3" sin a - 2{a'f'(3' cos a - {(3' f sin^ a cos a = 0, (4.81) 



I.e. 



P = ±[— p==+do (4.82) 

J sin a \/ d — sin a 



for some constants do ai^d (i > 1. Moreover, the system (j4.79p says that 

622 = ^(^ + |^)((a')^-(/3'sina)2), 

g M Q' , . 
023 = 7^( h -^)a 13 sma, 



(4.83) 



where 

Q = {a'f + {13' smaf . (4.84) 
Substituting (|i:83]) into diTTH]) . one finds that 5f(^ + ^) is a constant. We write 

5(^ + 1^) = 2c. (4.85) 

LP 2Q 

Thus 

ip = hQ^exp{J ^) (4.86) 
with /ci G M. Together with ()4.76p . one gets 

fi = hQ-'g-'eM[ (4.87) 
J 9 

with A;2 G M. Hence 

Q' c 



- 4g , 



a'a" + (3' sin a (/?' sin a)' c 
2((a')^ + (/3' sina)2) 5' 



(4.88) 



Theorem 4.4. Let a and g be arbitrary functions oft. The functions (3, Q, ip, a22, 033; 0^23? A*? 
ir, A, (ps andi^s are given by Km . (ITOI) . KW . KTSh . (jM) . (I07D . (ICT) . (IM) 

and (I4.57p . respectively. Then we have the following solutions of the MHD equations (ll.ip - (ll.4p ; 



w =(— 



(— — ) cos a + a' sma)X — (a' cos a + 022 sma)3^ 

4Q 5 

/?' cos a + 023)^ sin a — tpfiQCr'Ps cos a =F ^IJ'Q^(,r'4'sCt' sin a, 

-— ■ ) sin a cos /3 + a' cos a cos — 3' sin a sin /3) A" 

4Q g 

+ ( — a' sin a cos /? + 022 cos a cos /3 + 023 sin (3)y 
+ (— /?' sin^ Q cos /3 + 023 cos a cos /? — 033 sin P)Z 

ipi^iQCr't's COS /? sin a it {a' cos a cos/3 — /?' sin q sin (3)ipfiQ^ ^r'>Ps 
Q' c 

— — ) sin a sin /3 + a' cos a sin + sin q cos I3)X 

4Q 5 

+ { — a' sin a sin /? + 022 cos a sin /? + 023 cos I3)y 
+ [—(3' siv? a sin [3 + 023 cos a sin (3 — 033 cos 

— ipiJ.QCr't's sin a sin /3 it (a' cos a sin + sin a cos P)ipijLQ^ ^r'ips 

2c 4c 

— — sin a((a')^ — sin a)^)y — 7^0'/?' sin^ 

I 1 

— g(pfiQCr4>s cos a =F gft-i-d Q'^S.ri^s sin a, 



,2c 



4c 



H =(— cos a cos/3((a') — (/?'sina) ) — — a'/3' sin a sin/3)3^ 



Q 



4c 2c 
+ (— sin a cos a cos P + sin /3((a')^ — (/?' sin a)'^))Z 

— gipfiQCr<Ps sin a cos (3 ± {a' cos a cos j3 — (3' sin a sin (3) gipfiQ'2^j.ipg^ 

2c 4c 
i?^ =(— cos a sin /?((a')^ - (/?' sin a)^) + -pra'P' sin a cos /3)3^ 



Q 



Q 



4c 2c 
+ {—a'P' sin a cos q sin (3 — -pr cos /?((a')^ — (/?' sin a)^))Z 

— giffiQCr^Ps sin a sin /3 it (a' cos a sin (3 + 13' sin a cos (3)gipfiQ^(,r'ips, 



(4.89) 



(4.90) 



-P(.{a'n + ail - <9) — + («22 + ^22 + 43 " 2a23/3' cos a - (a'f)^ 



2 



(033 + 0^3 + a23 + 2023/?' cos a — (/?' sin a) 

— {a + 2a'a22 + 2023/?' sin a — (/3')^ sin a cos a)Xy 

-{{13' sin a)' + 2a'a23 + a' (3' cos a)XZ 

+(023 - aii023 + (022 - 033)/?' cos a - a'/3' sin a)3^2^ 

+{5r,i + 4a<5,,o)y'V'Q'^ + (<5.,i + 4M,,o)¥^V'Q'f 



2/i^r 



-L^ rF>(9' 2c Q' /— , a" /3'sina , 

±(-1) /iV<3(— H ^ TTTTiCrV's T <ffiy/Q{a { h — 7 + 022 — (023 - P cosa))3^ 

g g 2Q ip a' a' 



a 



ni ■ 1^' (/3'sina)' ^ , ^, , n 

+13 sn\a{ V — — h 033 + -— (023 + P cos q)Z)^,.V's), 

ip p' sm a p' sm a 



(4.91) 



where c is an arbitrary constant. 



Set 



5e7iy+722 _ ce-(7i3^+722) if r = 1, 
csin(7i3^ + + 6) if r = 



(4.92) 



, 5p7iy+722 _|_ „p-{7iy+722) :r 1 

C.= , , ' , .^ ' (4.93) 



and 




and 

. -/ 

ccos(7i3^ + 72 2^ + b) if r = 

where 71 and 72 are functions of t, and 6, c are real constants. 
Suppose that 

- (51 + (710- + 72T)Cr)\ 

qi + crir (4.94) 

Sl + T^r ) 

{92 + (TiC^I + 72Tl)Cr)\ 

92 + fTiCr , (4.95) 

S2 + n& ) 

where A = {aij)3x3 and B = (6ij)3x3 are 3x3 matrices whose entries are functions of t; the 
functions fj, qj, sj, a, T,ai and ri are also functions of t to be determined. 
The equations ()4.16p become 

b'liX + b[^y + b[^Z + - + (7ifTi + 72ri)';f )Cr 

- i-^Yiliy + 72^)(52 + ilKTi + 72n)^)^r + (a'3^ + Zp'sma){bn - ilicJi + 72Ti)Cr) 

+ ((021 - a')'^ + ^223^ + (023 + ^'cosq)Z + qi + o-^r)(6i2 - (-l)''7i(5'2 + (71 ^^l + 72n)'^)Cr) 

- ((/?'sina - 031)^" + (/3'cosa - 032)3^ - 033^ - si - r^r) 

X (^13 - (-l)"" 72(52 + {llCri + 72n)'V)Cr) 

+ ((^22 + 033) + (710- + 72'7")Cr)(bll'^ + bi2y + 613^ + /2 - (^2 + (7lf^l + 72n)'Y)Cr) 

- (a' + ai2 - (-l)''7l(s'l + (71 fJ + 72T)^)Cr)(&21'^ + ^223^ + &23^ + Q'2 + f^lCr) 

- (/3'sina + 013 - (-l)"" 72(5^1 + (710- + 72T)A')^r)(^31'^ + ^323^ + ^33-2^ + S2 + (J2ir) 

- (-fell + (7iO"i + 72Ti)Cr)(aii^ + ai23^ + ai3^ + /i - {qi + (71 + 72r)A')^,.) 

+ {-Ifilill + 7|)(52 + (7i'7i + 72^)^")^ = 0, (4.96) 



and 



b'^^X + 6'223^ + ^23^ + 92 + ^'iCr + (7l3^ + 72^)^^ier 

+ 621(011;!^ + (ai2 + a')y + (ai3 + /?' sin a)Z + /i - (51 + (71^7 + -f2T)X)ir) 

+ (&23 + 720-1 ^r) ((031 - P' sm.a)X + (032 - [3' COS a)y + 033^: + Si + T^r) 

+ [a' - a2i)(6ii;f + 6123^ + 613^ + /2 - {92 + (71^1 + 72n)'^)Cr) 

- (022 + ll(yir){b2lX + 6223^ + ^23^ + 92 + CFlir) + (^22 + ll(^lir){-Oi' X + Z0 COS O) 

- (023 + ^' cos a + 720-Cr)(fe31'^ + ^323^ + ^33^ + S2 + nCr) 

+ (&22 + 7if^ier)(a2i^ + 0223^ + 023^ + <?i + oi,) - (-l)'^r?(7? + ll)<y\ir = (4.97) 



bi,^X + 6^523^ + &33^ + S'2 + t{C. + (7^3^ + 72^)TlCr 

+ 631 {auX + (ai2 + a')y+ (013 + /?' sin a)^ + /i - (51 + (7if7 + j2T)X)Cr) 

+ (^32 + 7inCr)((a2i - a')X + a223^ + (a23 + (3' cos a)Z + qi + a^r) 

+ (^33 + 72TiCr)((a3i " sm.a)X + (032 - P' COS a)y + 033^: + si +r^r) 

+ {0 sina - a3i)(6ii;f + 6i23^ + 613^ + /2 - (52 + (71 '^i + 72Ti);i^)^r) 

+ {13' cos a - 032 - 7lT,^r,)(621'Y + &223^ + &23^ + Q'2 + (^lir) 

- (033 + 72rer)(fe31'^ + ^323^ + 633^ + S2 + n^^) - {-If v{ll + 7l)Tier = 0. 

(4.98) 



The coefficients of f^Cr and say that 

{Tia — TfJi = 0, , ^ 

(4.99) 
(710-1 + 72x1)51 = (71 cj + j2r)g2- 

Moreover, the coefficients of ^rCr in ^iz = ^3X yield 

7iT = 720-- (4.100) 

The above two equations suggest that 




where 



(4.101) 



h = —. (4.102) 

gi 



Motivated by Theorem 4.4, we assume that 

021 = —0-12 = ct', aai = — ai3 = /3'sina. (4.103) 
The coefficients of in ^ly = ^2X and the coefficients of y in (j4.96p show that 

bii = bi2 = bi3 = 0. (4.104) 
The coefficients of and Cr in (05]) - (l05D yield 

621 = 631 = 0. (4.105) 

Now 

Ri = (a'li + afi + (/?' sin a)'^)X — (a" + 2a'a22 — 2a32/3' sin a + (/3')^ sin a cos a)3^ 
+ (—(/?' sin a)' — a'{a23 + /?' cos a) — 2033/?' sin a — a'a23)Z 

+{-2a'a - 2/?'r sina - (-l)'"(5'i + (710- + 72r)^)[(7i + 71022 + 72(^32 - P' cosa))y 
+(72 + 71 (023 + (3' cos a) - 2033/3' sin a)Z + 71^1 + 72Si]}Cr + {-(^'i + (7io- + 72^)'^^) 
-an (51 + (710- + 72^")'^) - {aiiX + /i)(7ifT + 72T) 
+ (-l)^z.(7? + 72)(5i + (71^ + 72r)^)}Cr 
+{Abc5rfi + c^(5r,i)(7if^ + 72t)(s'i + (71 fJ + 72r)^) 

+f[ + aii/i - 2a'gi - 2siP' sin a, (4.106) 

i?2 = (a" + 2a'aii — (/?')^ sina COS a)A' 

+(022 + 022 + ^23(032 - cos a) - (a')^ - 032/?' cos 0)3^ 

+(^23 + 022(023 + /?'cosq) + 023033 — a'P' sina — 033/3'cosq)Z 

+ {a' + 022cr + (023 - /?' cos a)r - {-iy'i^{'yf + J2)^)^r 

+{-2a{gi + (71 (7 + 72t)^) + a[{'y'i + 71O22 + 72(032 - /?' cos a))y 

+ (72+71(023 +/?'cosa) - 2o33/?'sina)Z + 7igi + 72Si]}Cr 

+(t(7i(T + 72r),^rCr + q'i + 02291 + (023 - /?'cosq)si, (4.107) 

^3 = ((/?' sina)' + 2oii/3' sina + a'/3' cos a)^%' 

+(032 + 032O22 + 033(032 - /3'cosa) - a'/?' sin a - 022/?' cos a)>' 

+(o'33 + 032(023 + /3' cos a) + O33 - (/?' sin a)^ + 023/?' cos a)2^ 

+(t' + 033r + (032 + /?' cos a)a - {-lYiy{jl + ^2)^% 

+{-2{gi + (7i(T + 72T)A')/3'sina + t[(7i + 71O22 +72(032 - P' cos a))y 

+(72 + 71(023 + P'cosa) - 2033/?' sin a)Z + 71^1 + 72Si]}Cr 

+r(7i(T + 72'r).^,.Cr + s'l + 033S1 + (032 + /?' cos a)qi + 2/i/?' sin a, (4.108) 



Gi = (-ir/i'(7ifT + 72r)(9i + (7if^ + 72T);i^)C' 

+ {-lYh{gi + (7ifT + 72T)X){-fib23Z + 72^323^ + 7192 + 72S2)Cr-, (4.109) 



G2 = b32ib32-b23)y + h\-lYji{gi + {jia + j2r)XyCrC 



r 



+h{{b32 - b23)T - i-lYf27ligi + (7lf^ + l2T)X))Cr + (^'32 " 623)^2, (4.110) 

G3 = 623(632 -&23)^ + /i'(-l)'^72(5l + (7l^T + 72r);f)2erCr 

-h{{b32 - b23)(T - (-1)7271(51 + (71 + l2T)X))Cr. (4.111) 



By (j4.18p . we have that 



a" (/3')^ sin a cos a /3'smaa32 

033 = -, \ -, • 

a a a 

(/S'sina)' a' cos a a'a23 
022 = -T^r-. \ : h 



(5' sin a sin a f3' sin a 



and 



632 
623 



/^0^72 

72 + 71 (023 + 1^' cos g) + 72033 
/io/171 

But the coefficients of ^ in (j4.96p yield 

^(7i +71^22 + 72(032 -/?'cosa)) 



6: 



'32 



^23 



72 

^(72 + 7i(«23 + 13' COS a) + 72033) 



71 



Thus 

623 = 632 

and 



71 + 71^22 + 72(032 - /?' COS a) = 0, 

72 + 71(023 + (i' cos a) + 72033 = 0. 



Again by (j4.18p . we have that 

+ 2an-(-l)V(7i^ + 72^) =0, 



f (7if^ + 72T)' , , -,^r,,. 2 



71 <7 + 72T 

a'l , , /i(7iC^ + 72t) I ..r ( 2 , 2n , (-l)>o/2^(7l'7 + 72t) 

— + an + (-1) i/(7i + 72) + = 0, 

91 91 91 



(4.112) 



(023 - 032)011 = 023 - "32- (4.113) 

The coefficients of Cr in ^ly = ^2x say that 

7[ + 7ia22 + 72(032 - cos a) 



(4.114) 



(4.115) 

(4.116) 
(4.117) 



(4.118) 

7i/3' sin a = a'72 (4.119) 



and 

7i'?i + 72^1 + Ato^(7i92 + 72S2) = 0. (4.120) 

Since 

/i(7i9i + 72'Si) + 71 92 + 72S2 = 0, (4.121) 

comparing with (|4.96p . we get 

71 91 + 72S1 = 71 92 + 72S2 = 0. (4.122) 



Furthermore, by (j4.1UU|) and (|4.119|) . we can write 

1^^=^"; and 1" = ^"; (4.123) 

for some functions if and /i. 

From (lOGll - dOSi ). we have that 

f (/i(7i^ + 72t))' - (-l)^r?/i(7f + 7l)5i = 

\ {hgiY - haiigi + hfi{-yia + j2r) - /2(7i'7 + 72t) - 77/1(71 +72)51 = 0, 

J {ah)' - auha - hT{a23 + /?' cos a) - (-1)'' 77/1(71 + 72)0- = 
\ {rh)' + (/?' cos a - a32)ha - a^shr - (-1)'' 77/1(7^ + j^)t = 

and 

/2 - ail/2 = 



(3' sin a ip a /3'sina /3'sina (/J'sina)^ 



Thus 



The system ()4.118p says that 
where 

Then we get 



(4.124) 
(4.125) 



92 ~ o,22q2 - (023 + P' cos a)s2 = (4.126) 

^ §2 + iP' COS Q - a32)g'2 - 03352 = 0. 

For simpHcity, we assume that 

f2 = q2 = S2 = 0. (4.127) 

Then by (|4.112p and (|4.117p . one has 

(P' sin a)' a' 13' cos a a'a23 

022 = — \- — h 



(5' sin a /?' sin a /?' sin a 

(/?'sina)' a'/3'cosa a'a23 -,r,o^ 

^33 = 75r^ (4.128J 

(/9 p' sm a p' sm a p' sm a 

a' ,(p' a" (/3'sina)' a'/?' cos a, ^, (q;')^023 
az2 = --^. {— + — + ^^- - + ^r-- ) + /3cosa 



an = — . (4.129) 



3^ + ^ + % = i-iyi^^Q, (4.130) 
ip fj, Q 

Q = {a'f + {P' sin af. (4.131) 



^ = exp( y (-l)V(^Qdt) (4.132) 



for some constants Ci. 

We consider the special case of (j4.128p . in which 

023 = 032- (4.133) 

Hence we have 

a'B'sina,(p' a" (B'sina)' a' 3' cos a (/3')^ sin a cos q, 

023 = ^ - + — + . ' + -ir- — 1 • 4.134 

Q if a p' sm a p' sm a a' 

One shows that the system (I4.124p is compatible by ()4.118p . Moreover, the equation (l4.125P i is 
compatible with (j4.124p if and only if 

a" (3' sin a - a' 13" sin a - 2{a'fl3' cos a - {I3'f sin^ a cos a = 0, (4.135) 



I.e. 

P = ± [ — , _^ + do (4.136) 

J sin a y d — sin a 



for some constants do and d > 1. In this case, one shows the system (I4.125P is compatible. Thus 
one gets that 

h = C2(^^ exp((-l)'^ y"(r/ - z/)(^Qdt), (4.137) 



where C2 is a real constant. 
Write 

I s\ = \a' , 



qi = X(3' sin a. 



(4.138) 



Then we have 



+ (— - Q)^ + (022 + 022+0^3 - 2a23/?'cOSQ - {o!f) — 

+ (o'as + 023 + O33 + 2a23/3' cos a — {j3' sina)^)— + (a" + 2a' (/3')^ sin q cos a)Xy 



2 ' (f 

+ ((/3'sina)' + 2^/3'sina + a'/3' cos a)XZ 



+ (023 - 023— + (022 - 033)/?' cos Q - a j3' cos a)yZ 

+ -{Ahcdrfl + c25,,i)(5i + (/j^Q;I^)2 + (/( + ^ - 4Aa'/3'sina);i^ 
^ Ji 

+ ((A/3 sin a)' + 022 A/3 sin a + (023 — /3' cos a)(Aa'))3^ 

+ ((Aa')' + 2/i/3' sin a + (023 + /3' cos a) A/3' sin a + 033 Aa')2^ 

- 2^xQXir-2^ + {-lY^{^-^-{-lYv^^Q)Cr) (4.139) 
modulo the transformation in (jl.9p . where C2 is an arbitrary constant. 

Theorem 4.5. iygi a, gi, \, fi be functions oft. The functions fj,, Cr, P, Q, 0-22, 033, 023 
are given by KT32\i . K92\i . KU3\\ . lKTm\\ . IKTEB . (ITOSD and (gJMD . respectively. Then we 
have the following solutions of the MHD equations (jl.ip - (jl.4p .' 

u =( — cos a — a' sin a) A' — (a' cos a + 022 sinQ)3^ — (/3' cos a — 023)-^ sin a 

+ /i cos a — A/3' sin^ a — {gi + (ffiQX)(r cos a — fi^rOt sin a, 

=(— sin a cos /3 + q' cos a cos fi — 13' sin a sin /3).%' 

+ (—a' sin a cos /3 + 022 cos a cos /3 — 023 sin j3)y 
+ {—13' sin^ a cos /3 + 023 cos a cos /3 — 033 sin /3)^ 

+ /i sin a cos [3 + X(3' sin a cos a cos /3 — Aa' sin /3 (4. 140) 

— {gi + ippQX)C,r sin a cos /3 + ^a'^,. cos a cos /3 — /x/?'^r sin a sin /3, 

w =[— sin a sin /3 + a' cos a sin (3 + (3' sin a cos /3) A' 

+ (—a' sin a sin /3 + 022 cos a sin [3 + 023 cos /3)3^ 
+ (— /3' sin^ a sin [3 + 023 cos a sin /3 + 033 cos I3)Z 
+ /i sin a sin /3 + A/3' sin a cos a sin /3 + Aa' cos (3 

— {gi + ippQX)C,j. sin a sin j3 + /ia'^^ cos a sin /3 + /x/?'^,. sin a cos /3, 



= — C2(/3^ exp((— l)*" y ~ '^)'-pQ'ii)i{gi + fl-''Q'^)Cr cos a + i-La'^r sin a), 

=C29?^ exp((— 1)*^' J ~ '^)fQ'^t){—{gi + LpfiQX)(r sin a cos /? 

+ ^a'^r cos Q cos /? — sin a sin /?) , (4.141) 
=C2'P^ exp((— 1)'' y ~ i^)'pQ(it){—{gi + (pfiQX)Cr sinasin/3 

+ fia'S^r COS a sin /3 + /i/3'.^r sin a cos /?) 



and p is given by ()4.139p . 



Set 



e"i^ - ae-"i^ if r 



sin(ai3^) 



= 0, 
if r = 1, 



and 



be^^y - ce-^^y if r = 0, 

csin(ai3^ + 6) if r = 1, 

6ie"i^ - cie""!^ ifr = 0, 

ci sin(ai3^ + 6i) ifr = l, 



Cr-- 



oaiy 



+ ae-^^y if r = 0, 



cos(ai3^) 



if r = 1, 



6e"i^ + ce-"i^ if r = 0, 
ccos(ai3^ + 6) if r = 1 



ci cos(ai3^ + 6i) ifr = l, 



where oi is a functions of t and a, b, bi, c, ci are real constants . 
We assume that 

U = auX + ai2y + aisZ + aif{6r + a2y4>r - ^Cr - aia2X'^ipr), 

V = a2\X + a22y + 023^^ + !{ir + 1axa2X(^r), 
W = a3i;f + 0323^ + 033^ 



and 



V} = hiiX + bi2y + 613Z + aig{er + 623^0r - XCr - aib2X^i;r), 



V = b2lX + b22y + 623^ + giCr + 2aib2X<Pr), 
W = bsiX + bs2y + ^33^, 

where aij, bij, a2(/ 0), 62(7^ 0), /(/ 0) and g{^ 0) are functions of t, and 

J an + 022 + 033 = 0, 
\ ^11 + &22 + &33 = 0. 

Substituting (I4.145|) - (I4.147I) into KWh and KT8\i . we get that 



(4.142) 
(4.143) 

(4.144) 



(4.145) 



(4.146) 



(4.147) 



Ri = a[^X + a'^^y + a[^Z + (ai/)'(0,. + a23^0r - ^Cr - aia2;fVr) 

+ «i/(a'i3^er + a23^<Ar + a[a2y^iJr - i-iya[Xy^r - (aia2)''^V'r - {-!)'' a[aia2X^y(j)r) 

+ (an - ai/(Cr + 2aia2X'tpr))ianX + (ai2 + a )3^ + («13 + P' sma)Z 

+ ai/(6'r + a23^</'r - '^Cr - aia2'^^V'r)) 

+ (ai2 + ai/(aier. + a2<l)r + aia23^V'r - (-l)''oi'^^r - (-l)''«ia2'Y^0r)) 

X ((a2i - a)X + a22y + (^23 + P' cos a)Z + /(^r + 2aia2'^'</'r)) 

+ «i3((o3i - /3' sin a)^ + (032 - /?' cos a)3^ + 033^^) 

- a'{a2iX + 0223^ + ^23^ + fiCr + 2aia2X(j)r)) - {a-iiX + 0323^ + 033^)/?' sin a 



i-iyualfiOr + a2y(pr - XQr - aia2;fVr 



(4.148) 



Gi = {b2iX + 6223^ + h2zZ + g{ir + 2a^b2X(l)r)) 
X (^21 - bi2 + aig{-aier + b2(t>r - aib2yil^r + (-l)^ai('^Cr + aib2X'^ (t>r)) 
- {bn-bu){bziX + b^2y + b^^Z), (4.149) 



R2 = a2iX + a'22y + + f'iir + 2aia2X4>r) 

+ f{a[yCr + 2{aia2yXcj)r + 2a[aia2Xy^r) 

+ (021 + 2aia2f(l)r){auX + {au + a')y + (013 + /3'sma)Z 

+ aif{9r + a2y4>r - XCr - aia2'^^^r)) 

+ (a22 + aif{Cr + 2ala2X^lJr)){ia'2l - ol)X + 022!)^ + (023 + cos a)Z + /(^^ + 2aia2X(pr)) 

+ a23((a3i - /3' sin a)X + (032 - /3' cos a)y + 033^) 

+ a'{auX + ai2!y + 013^ + aif{9r + a2y(t>r - X(r - 0102^'^^/;^)) 

- [asiX + 0323^ + a3sZ)P' cos a - {-ly'ualfi^r + 2aia2X^r), (4.150) 

G2 = (632 - b23){b3iX + 6323^ + 633^) 

- (^21 - 612 + aigi-aiEr + b2(f>r - ai62>'^/'r- + (-l)''ai'^Cr + (-l)''ai^2'^^<^r)) 

X {bnX + bi2y + bi3Z + aig{9r + b2y(l)r-XCr-aib2X^i;r)), (4.151) 

i?3 = ag;^ A" + a32>' + 033^ + 031 {aiiX + {au + a')y + (ai3 + /3' sin q)Z 
+ aifiOr + a23^0r - XCr " aia2^Vr)) 

+ 032 ((a2i - a')X + a223^ + (023 + P' cos a)Z + /(^^ + 2aia2X(pr)) 

+ 033 ((^31 - P'sma)X + (032 - /3'cosa)3^ + 033^^) 

+ {aiiX + 0123^ + ai^Z + aif{9r + a23^<Ar - '^Cr - aia2X'^\(jr))P' sina 

+ {aiiX + ai2y + anZ + f{Cr + 2aia2X (l)r))P' cos a, (4.152) 

G3 = (fel3 - ?>3l)(&ll^ + ^123^ + 613^ + + ^2J0r - ^Cr " 0162^-^^)) 

- (632 -623)(?>2i^ + 6223^ + 623^ + <7(er + 2ai62^</'r.)), (4.153) 

ft'llA" + 6;23^ + b'lS^ + iaigyi0r + &23^<Ar " XCr " aifo2^Vr) 

+ ai(7(a'i3^e, + 6^3^</>, + a;623^Vr- - (-l)"a'iA'3^er- - (ai62)''^V'r - (-l)"a'iai62^^3^</'r) 

+ {a'y + Z (3' sin a){bu - aigiCr + 2alb2X^Pr)) 

+ ((^21 - a')^ + a223^ + (023 + P' cos a)Z + /(^^ + 2aia2X^r)) 

X (612 + aig{aier + b2(t>r + aib2yipr - {-'^Yai{XS,r + aib2X'^ (pr))) 

+ ^i3((«3i - l3' sin a)X + (032 - P' cos a)y + 023^^) 

+ {a22 + a33 + f{aiCr + 2ala2Xtpr)) 

X (6ii;f + fei23^ + bisZ + aic/(er + b2y^r - XCr - aib2X'^^r)) 

- {a + ai2 + aif{aier + a24>r + aia2yipr - {-'^Yai{XS,r + aia2X'^ (pr))) 

X (621^ + 6223^ + 623^ + <7(er + 2ai62^(/'r)) - (013 + /3'sina)(63i^ + 6323^ + 633^) 

- {auX + 0123^ + anZ + aif{9r + a23^0r- - XCr - aia2X'^ipr)) 
X {b22 + bss + g{aiCr + 2alb2X'il!r)) 

- i-lYr]alg{9r + b2y^r - XCr - ai62A'Vr) = 0, (4.154) 



b'^^X + b'^^y + 6'23^ + g'i^r + 2aib2X(t)r) + ^(a'lXr + 2(ai62)''^'<^r + 2a\aib2Xi;r) 
+ (621 + 2aib2g4>r){aiiX + (ai2 + a')3^ + (013 + (3' s\na)Z 

+ ai/(0,. + a2y<t)r - XCr - aia2X'^ijr)) + {-a'X + cos a)(622 + aig{Cr + 2ai62'^'V'r)) 
+ &23((a3i - (3' suia)X + (032 - 13' cos a)y + 033^) 

+ (a' - 021 - 2aia2f^r){biiX + bi2y + 613^ + ai5(^r + b2y^r - XQ- - aib2X^i:r)) 

- (a22 + ai/(Cr + 2aia2X^r)){b2iX + 622^^ + ^23^ + g{& + 2aib2X(^r)) 

- (a23 +/3'cosa)(63i'^' + fe323^ + ^33^) 

+ (622 + aig{Cr + 2aib2X%br)){a2iX + a223^ + 023^ + /(^r + 2aib2X^r)) 

- {-lyrjalgi^r + 2aib2X(f)r) = (4.155) 

and 

b'si^ + b'32y + b'ssZ + &31 {anX + {a^ + a')y + (023 + /3' sin a)Z 

+ aif{9r + a2y4>r - X(r - 0102^"^?/;^)) 

+ &32((a2i - a')X + 0223^ + (a23 + P' COS a)Z + /(^^ + 2aib2X^r)) 

+ (/3'sina - a3i)(6ii;f + buy + bi^Z + ai5(0,. + 623^</'r - -^Cr - aib2X^^r)) 

+ (/?' COS a - a32)(62i'^ + ^'223^ + ^23^ + g{ir + 2aib2X cpr)) " a^^ib^iX + 6323^ + 633^) 

+ 633(a3i;i^ + a323^ + a33^) -^'33('^'/9'sina + 3^/3'cosa) =0. (4.156) 

From the coefficients of A'^c^j. in (^iz = $3A'i we get 

(^23 + cos 01)02/ + /U0623&25 = 0. (4.157) 

While, from the coefficients of X'^cjiy in ^'i = 0, we have 

- (a23 + 13' cos 0)625 + 02/^23 = 0. (4.158) 

Since the functions 02, 62, / and g are nonzero, we obtain that 

{ 023 = — cos a, 
623 = 0. 

The coefficients of 4>r in = suggest that 

032 = -/?'cosa. (4.160) 
The coefficients of C,r in ^2Z = '^3:>' say that 

632 = 0. (4.161) 

Thus the coefficients of in ^3 = imply 

(3' = 0. (4.162) 

Substituting (|4.159p - (|4.162p into (|4.148p - (|4.156p . we can simplify the calculation. One shows 
that 

{ai3 = 031 = 0, ^ ^ 

(4.163) 
bij = for I, j = 1, 2, 3, 

I = = = ^' (4.164) 

[ an = ^22 = ^33 = 0, 

and 

i-^Ymlt (4.165) 



(4.159) 



g = C2e 
62 = 02 = constant, 



where fc, ci and C2 are constants. Thus one easily shows that 

p = p{_a2e(-i)'2i'»?t[l(45jQ^^ + c5i,r)cX'^ + a2i26o,r{ab + c) + 6i,rCcosb)X^ 
+ a'la2{4:a5Q^r + Si,r)-^ — h a2(2(5o,r(a^ ~ c) + 6i^rCsinb)Xy + (5i,rCCi sin(6 — bi) 
+ 2(5o,r(bci - bic))X + ^(2(5o,r(ci - a6i) - (5i,^ci sinbi)^'^] + —Af^ 

X a2e(-l)''^'^?*(aia23^(^r^r- - y (ai(4a(^o,r + 5l,r)3^^ + ^r) + ^rCr 



(<5o,r2(a6i + ci) + 5i,.ci)3^) + A;(e, + a23^^. - ^</..)] - ^e'--'r2u4t^2 ^ ^y2^ 

ai 2 2 

/io/o{aiC2e("^)'^''''i*(-ai^rer + a2^r<Pr - aia23^^r^r)'^ 



+ ((— l)''ai^^ - 2aia2<?!)rer + 2aial(j)f. - 2alaly<j)ripr)-^ 

+ i-iyaia2X%<Pr + i-iyalal^^l - '^C2e^-^y^^-l\er + a2y<Pr? - 2a\a2C2e^-^^'''^''-' 

X ('^o,r[77^(&ie2"i^ - - (ci + a6i)3^] + sinfti sin2(ai3^) 

, ci ai3^ sin(2ai3^) 
+ -cos6i(— + ^ )]) 

- 2alc2e^-'^'''^^'\6o,r[-aly^ + ^(e^-^^ - e'^-^^) - ^(e^'^^^ + e'^'^^^)] 

. r5 sin2(aiy) aiysm^jaiy) aiy 

- ^i,r[g sin(2ai3^) + — ])} (4.166) 

modulo the transformation in p.9p . 

Theorem 4.6. Let k, I, ai, a2, Ci and C2 be any real constants. Then we have the following 
solutions of the MHD equations (jl.l|) - p.4p .- 

' u =kX sm{kt + l)-ky cos{kt + I) - cie(-^)''^"?*(^^ + 2aia2Xcl)r) sm{kt + /) 

+ aicie("^)''''"i*(6lr. + a2y(pr - ^Cr - aia2X^Tpr) cos{kt + I), 

V =kXcos{kt + 1)- kysm{kt + I) + aicie("^)'''"'?*(6',. + 023^0^ - ^Cr - aia2X^il^r) 

X sin(A;t + /) + cie^''^^'' ""^^^^r + 2aia2X(j)r) cos{kt + /) 
I, w =0, 

(4.167) 

' =aiC2e(-^)'''''i*(^r + a2y(l)r - ^Cr - aia2X^^r) cos{kt + 
- C2e(-^)''''''i*(Cr + 2aia2X(pr) sm{kt + I), 
if2 =aiC2e(-^)'''''i*(0r + 023^0,- - XCr - aia2X'^i;r) sm{kt + /) (4.168) 

+ C2e(-^)''''''i*(Cr + 2aia2X(j)r) cos{kt + I), 
H'^ =0 

and p is given by (j4.166p . 
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